Nonequilibrium wetting transitions with short range forces 
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We analyze within mean-field theory as well as numerically a KPZ equation that describes nonequi- 
librium wetting. Both complete and critical wettitng transitions were found and characterized in 
detail. For one-dimensional substrates the critical wetting temperature is depressed by fluctuations. 
In addition, we have investigated a region in the space of parameters (temperature and chemical 
potential) where the wet and nonwet phases coexist. Finite-size scaling analysis of the interfacial 
detaching times indicates that the finite coexistence region survives in the thermodynamic limit. 
Within this region we have observed (stable or very long-lived) structures related to spatio-temporal 
intermittency in other systems. In the interfacial representation these structures exhibit perfect tri- 
angular (pyramidal) patterns in one (two dimensions) , that are characterized by their slope and size 
distribution. 

PACS numbers: 5.10.-a,64.60.-i,68.08.Bc 
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I. INTRODUCTION 

When a bulk phase a is placed into contact with a sub- 
strate, a layer of a second, coexisting, phase (3 may form if 
the substrate preferentially adsorbs it. At a wetting tran- 
sition, the thickness of the (3 layer diverges. Equilibrium 
wetting has been experimentally observed and theoreti- 
cally investigated using, amongrnany other techniques, 
interface displacement models Qj 0, Q • Within this ap- 
proach one considers the local height of the a(3 interface 
measured from the substrate, h(x), and constructs an ef- 
fective interface Hamiltonian, H.(h) 4]. In equilibrium 
situations, one typically has 



H(h) 



v{Vh) 2 + V{h) 



(1) 



where v is the interfacial tension of the af3 interface (or 
the interfacial stiffness if the medium is anisotropic) and 
V(h) accounts for the interaction between the substrate 
and the a(3 interface. 

If all the microscopic interactions are short-ranged, one 
may take for sufficiently large h at bulk coexistence Q 



V(h) = b(T)e- h + ce" 



2h 



(2) 



where T is the temperature, b(T) vanishes linearly as 
T — T w , T w being the wetting temperature, and c > 
5]. By minimizing Q one finds Q Q a critical wet- 
ting transition at b = 0, i.e. the interface height (or 
equivalently, the wetting layer thickness), (h), diverges 
continuously as b ~ T — Tw — * bw = — . Equilibrium 
critical wetting has been studied for decades and a rich 
(non-classical) behavior predicted Q, 0] . 

Wetting transitions may also be driven by the chemical 
potential difference between the f3 and a phases, fi. In 
this case wetting occurs at any temperature above Tw 
(i.e. for b > bw) as /z = is approached from the a phase. 



This is always a continuous transition and it is known 
as complete wetting 0. A study of complete wetting 
requires adding a linear term [ih to the Hamiltonian 

A dynamic model for the growth of wetting layers has 
been proposed through the Langevin equation 
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dt h(x,t) = - S -^-+r ! = isV 2 h-^-+r l , (3) 
where n is Gaussian white noise with mean and variance 



(Tj(x,t)> = 0, 

(77(x,i)7/(x',i')) = 2D5(t - t')S(x - x'). 



(4) 



Equation (|3 is an Edwards- Wilkinson (EW) growth 
equation yjHj in the presence of an effective interface po- 
tential V(h). It describes the relaxation of the interfacial 
height h towards its equilibrium value, i.e. the value of h 
that minimizes TL. Within this context \x can be viewed 
as an external driving force acting on the interface. Re- 
call that in the absence of the wall, the corresponding 
equilibrium states for [i < and fi > are the a and (3 
phases, respectively, whereas phase coexistence occurs at 
(i = 0. 

Equilibrium models, however, are not sufficient to 
study wetting in Nature, since in a wide range of phe- 
nomena (e.g., growth of thin films, spreading of liquids 
on solids or crystal growth) thermal equilibrium may not 
hold. Nonequilibrium wetting transitions have been re- 
cently studied by Hinrichsen et al. 9] in a lattice (re- 
stricted solid-on-solid) model with dynamics that do not 
obey detailed balance. The continuum nonequilibrium 
counterpart of this discrete model, is a Kardar-Parisi- 
Zhang (KPZ) equation in the presence of a bounding po- 
tential, whose properties have been anal yze d extensively 
by one of the authors and collaborators [ll], [l3j . Clearly 
this is the most natural extension of the EW equilibrium 
growth model to non-equilibrium situations. In fact, in 
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the absence of a substrate the KPZ non-linearity, \(\7h) 2 , 
is generically the most relevant noncquilibrium pertur- 
bation to the equilibrium EW equation The KPZ 
nonlinearity is related to lateral growth, and although 
this mechanism is unlikely to be relevant in simple fluids, 
it may determine the wetting behavior of systems with 
anisotropic interactions for which the grow th of tilted in- 
terfaces depends on their orientation 14]. For instance, 
it has been shown that crystal growth from atomic beams 
is described by the KPZ equation [l^. From a theoret- 
ical point of view a key and ambitious task is that of 
developing general criteria to establish whether the KPZ 
non-linear term should be included in a given interfacial 
model. 

Related works published recently by Miiller et al. [lfij . 
Giada and Marsili ^tJi Hinrichsen et al. and our- 

selves [l0| , consider similar non-equilibrium models in the 
presence of various types of walls. 

In this paper, we further study the KPZ interfacial 
equation in the presence of different types of potentials, 
attractive and repulsive. We focus on the connection 
of the associated phenomenology with non-equilibrium 
wetting and depinning transitions. In particular, we 
will stress that the transitions called "first-order non- 
equilibrium wetting" in [l8j are not wetting transitions 
but, rather, depinning transitions. Also, we study for the 
first time the two-dimensional version of this model, and 
report on new phenomenology. 

The remainder of this paper is organized as follows. 
In the next section we introduce our KPZ-like model. A 
mean-field picture is provided in section III, and its pre- 
dictions numerically tested in section IV. The conclusions 
are summarized in the final section. 



II. THE MODEL 

The model under study is defined by the Langevin 
equation 

d t h{x 7 t) = vV 2 h + \(Vh) 2 ~?X^+V, (5) 

oh 

where V(h) = -(a + l)h + be~ h + ce~ 2h /2, c > 0, rj obeys 
Q and a + 1 = /i is a chemical potential. 

In the absence of the exponential terms (the limiting 
wall) the interface moves with a nonzero positive (nega- 
tive) mean velocity for fi larger (smaller) than a certain 
critical value fi c . In one dimension fi c can be found an- 
alytically since both the KPZ and EW equations have 
the same Gaussian steady-state height distribution [8j- 
Thus /i c — X((\7h) 2 ), which for discrete lattices can be 
approximated by [n| He = — (D\)/(2vA), A being the 
lattice cutoff. Note that for A ^ a nonzero chemical 
potential is required to balance the force exerted by the 
nonlinear term on the tilted interfacial regions. For A = 
the model reduces to the equilibrium one and [i c — as 
usual. For negative values of A the interface is (on av- 
erage) pushed against the wall, while for positive A it is 



pulled away from the wall. Thus the behavior of the sys- 
tem is determined by the sign of A |20|. In this paper 
we will consider A < only (which corresponds to the 
case studied using microscopic models 0])- Results for 
positive values of A will be published elsewhere. 

It is our purpose to study the effects of a substrate 
that adsorbs preferentially one of the two phases on the 
stationary properties of the interface. This is achieved 
by considering b to be negative in equation ©. 

This same equation, J5J, has been recently studied by 
Giada et al. [l7j as a generic nonequilibrium continuum 
model for interfacial growth. However, their choice of 
control parameters [l7( privileges the role of the noise as 
the driving force of the nonequilibrium transitions. By 
contrast, motivated by the role of the chemical poten- 
tial and temperature in equilibrium wetting, we fix the 
noise intensity and choose a and b as control parameters 
that are the fields driving critical and complete wetting 
transitions. 

To establish the analogy with equilibrium, let us stress 
that just like in equilibrium complete wetting, nonequi- 
librium complete wetting occurs when the attractive po- 
tential V(h) is not capable of binding the interface, at 
temperatures above the wetting transition temperature, 
b > bw, as the chemical potential approaches that of 
'bulk coexistence', fi — > fi c . At this transition the in- 
terface begins to move and (h) diverges. On the other 
hand, the nonequilibrium analogue of critical wetting cor- 
responds to the unbinding of a bound interface at 'bulk 
coexistence', /x = fi c , as b — > b^y. 

In order to analyze equation JSJ it is convenient to per- 
form a Cole-Hopf change of variable /i(x, t) = — In n(x, t), 
leading to 

„ ^2 dV(n) . 

o t n = vy n h nrj, (6) 

on 

with V(n) = an 2 /2 + fen 3 / 3 + cn 4 /4. This describes the 
interface problem as a diffusion-like equation with mul- 
tiplicative noise [Til l2lj ]. In this representation, the un- 
binding from the wall ((h) — » 00) corresponds to a transi- 
tion into an absorbing state (n) — > 0. In the following we 
will use both languages, h and n, indistinctly although 
the natural description of wetting is in terms of h. The 
case b > was studied injft[ll|, while b < is the case 
studied in [H E IH 13 • 

Note that we have made use of Ito calculus, and 
thus equation (jSJ should be interpreted in the Ito sense 
[2^]. In general, potentials of the form bn p+2 /(p + 2) + 
cn 2p+2 /(2p + 2) with p > result in equivalent effective 
Hamiltonians since, when expressed in terms of h, p can 
be eliminated by redefining the height scale. The case 
p = 2 (with fixed b < 0) has been studied in (2^ in 
the context of stochastic spatio-temporal intermittency 
(STI). The unsuspected connections between these two 
problems are illustrated in the following sections. 
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III. MEAN FIELD 

In this section we analyze equation (jSJ) at the mean 
field level. We begin by discretizing ijfjj on a regular d- 
dimensional lattice 



d t rii 



v 
2rf 



dV{ ni 
dni 



rkVi 



(7) 



where rij = n(xi,t) and the sum over j runs over the 
nearest neighbors of i. The Fokker-Planck equation for 
the one-site stationary probability P(rii) can be easily 
worked-out. In mean-field approximation (i.e. substitut- 
ing the values of the nearest neighbors by the average n 
value), the stationary Fokker-Planck equation for D = 1 



is 



d_ 

On 



dV{ri) 
dn 



+v(n— (n))Pt(n) 



|l(n 2 P t (n))=0 (8) 



and its associated solution 



P(n, <n» = N\ ex P - f V'W + (n)) 

where the integration constant N is determined by a nor- 
malization condition and (n) is obtained from the self- 
consistency requirement 

In) = £± nP ^ <"» = F((n)) (10) 
{) I °°dnP(n,(n)) H ™ h {W) 

Let us consider two limiting cases where analytic so- 
lutions of Ea. (|10|l can be worked out. In the zero- 
dimensional case, or equivalently v — 0, the solution of 
© reads 



P(n) = N- 



exp - 



2p 



,a+2 



(11) 



that, in terms of heights, is P(h) ~ exp[(a + l)h — 
be~ h /p — ce~ 2h /2p] and yields the effective potential 



V e ff(h) = - \nP{h) = -{a + l)h + be- h /p + 



ce 



-2h 



Clearly this coincides with the potential in J5j|. A com- 
plete wetting transition occurs when approaching a = — 1 
with b > and critical wetting is found at a = — 1 with 
6 = 0. 

For spatially extended systems, in the v = oo limit, 
a saddle-point expansion in v yields V'(n) = (I?! l24| . 
Thus, the dynamical behavior in this limit is that of the 
deterministic mean- field version of (JHJ: for any p > 0, 
there is a line of second order wetting transitions at a = 
and b > 0, and a line of first-order transitions at a > 
and b = ~(p + 2)^/ac/(p+ 1). These lines meet at a 
tricritical point at the origin. 

For values of v other than zero or oo, the self- 
consistency equation (n) = F((n)) has to be solved nu- 
merically. Without loss of generality, we set p = 2, c= 1, 
and illustrate in Fig. nt ne three different regimes: one 



stable solution at (n) — (dash-dotted line); one unsta- 
ble solution at (n) = and a stable one at (n) ^ (solid 
line); two stable solutions and an unstable one (dashed 
line). Stable solutions can be identified by a negative 
slope of F({n)) — {n) at the intersection point [25|. A 
nonzero solution emerging continuously from (n) = as 
a function of a and 6, signals a second order transition. 
This is the case for the dash-dotted and the solid lines 
in the inset of Fig. ^ By contrast, when the nontrivial 
solution appears discontinuously as a function of a and 
b, the transition is first-order (dash-dotted and dashed 
lines). The corresponding phase diagram is depicted in 
Fig. ^ The solid line is a second-order phase boundary 
from non-wet to wet substrates. Between the two dotted 
lines, the wet and non-wet phases coexist as stationary 
solutions of the dynamical equation. The three lines join 
at the tricritical point at a = b = 0. 

In order to determine the order parameter critical ex- 
ponent in mean field approximation, we proceed as in 
[l3| . First, we rewrite (|rU|) as 

r°° / h t 2 P\ 

(n)- 1 = -d {n) In J dtt a cxp {-- tP - -^) e ' {n)t - 

(12) 

Next, we introduce a Gaussian transformation and ex- 
pand the resulting integrals for small (n). We find 
(n) ~ |a| 1 / p , and thus a c = and (3 = l/p. 



IV. BEYOND MEAN-FIELD THEORY 

In this section we explore whether the mean-field phase 
diagram structure survives when the effects of fluctua- 
tions are taken into account. Mean- field exponents are 
expected to hold above the upper critical dimension d c , 
which in the present case, equation ©, is known to be 
d c = 2, (corresponding to 3 bulk dimensions and in the 
weak coupling regime of the KPZ) For positive val- 
ues of b and d > 2, the second term n 2p+2 in the effective 
potential is irrelevant 01 an d the we are left with 



dtii = vV 2 n — an — bn p+1 



nr), 



(13) 



defining the multiplicative-noise (MN) universality class 

We have solved © numerically for different system di- 
mensions. In particular, for a one-dimensional substrate 
we have considered a system size L — 1000, v = p = 2, 
D = 1 and c = 1.5. The time step and the mesh size were 
set to 0.001 and 1, respectively. We started by determin- 
ing the chemical potential for which the free interface has 
zero average velocity. For the parameters given above we 
found a c ~ —0.064. Then we fix a — a c and calculate 
(n(t)) for different values of b and large t. Length and 
time units are in lattice spacings and Monte Carlo steps, 
respectively. 
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A. Critical wetting 

To study critical wetting we set a = a c and consider 
small values of b for which an initially pinned interface re- 
mains pinned, and increase progressively b until the non- 
wet phase becomes unstable at by/ . The critical point is 
estimated as the value bw that maximizes the linear cor- 
relation coefficient of log(n) versus log \b — bw\', the criti- 
cal exponent is then determined from the corresponding 
slope (Fig. |2J). It is found that the critical "tempera- 
ture" is depressed from its mean-held value bw — to 
b\y = —0.70 ± 0.01, with an associated critical expo- 
nent f3 — 1.20 ± 0.01 (the error in the exponent comes 
from a least-squares fit). Below (above) that value we 
find first (second) order depinning transitions, by vary- 
ing a. Therefore, as in mean field, there is a "tricritical" 
point, joining a line of second order transitions (b > bw) 
with one of first order transitions (b < bw)- The criti- 
cal exponents and universality of this multiplicative-noise 
tricritical point has not been investigated before. 

The finite coexistence region allows us to define critical 
wetting along a range of different paths, delimited by the 
dashed lines in the mean-field diagram of Fig. We have 
checked that the value of (3 does not change when the 
critical point is approached along different paths within 
this region. Further numerical and analytical studies of 
this new universality class will be left for future work. 



B. Complete Wetting, b > bw case 

We consider a one-dimensional substrate, with b = 1 > 
bw, let the system evolve to the stationary state and then 
compute the order parameter (n) for different values of a 
near its critical value. As a — > a c a continuous transition 
into an absorbing state n = is observed; it is the non- 
equilibrium counterpart of complete-wetting. The asso- 
ciated critical exponent is found to be (3 — 1.65 ± 0.05 
in good agreement with the prediction for the MN class, 
(3 = 1.5 ±0.1 Other positive values of b yield similar 
results. In addition, we have simulated systems above 
the upper critical dimension, in d = 3, with L = 25, 
6 = 5 and other parameters as in the one-dimensional 
case. Our best estimate for (3 is (3 — 0.96 ± 0.05, indi- 
cating that this transition is governed by the weak noise 
fixed point of the MN class [llj . For larger values of the 
noise amplitude we find a strong coupling transition, in 
agreement with the theoretical predictions 0, . Fi- 
nally, we note that both numerical and renormalization 
group arguments lead to = 1 in the weak coupling 
regime, independently of the value of p [llj. This is at 
odds with the mean-field prediction [3 = I /p. Birner et 
al. [? ] have recently suggested a transition from 1/p to a 
nonuniversal behavior depending on the ratio of the noise 
to the strength of the spatial coupling. However, this 
discrepancy appears to be generic since different types of 
mean-field approaches yield the same (incorrect) result 
and the origin of the discrepancy remains unclear |l3j . 



Since the Cole-Hopf transform of the MN equation is 
the same as KPZ with an additional exponential term, 
the MN exponents are those of KPZ iff the extra term 
is an irrelevant term of the KPZ renormalization group 
flow. Note that the Cole-Hopf transform fixes the value 
of \jv and thus the A = limit can not be considered 
when this transformation is used. In addition the poten- 
tial V(h) is a relevant term of the EW equation. In this 
regime adding a non-linear potential is a relevant pertur- 
bation and it does indeed change/determine the wetting 
exponents (cf. with the literature on equilibium wetting 
lj). Thus EW plus a (non-linear) wetting potential is 
not equivalent to KPZ in the weak coupling regime plus 
the same wetting potential. 



C. Depinning transition at b < bw 

As expected, no transition is found as a — > a c when 
equation JBJ is solved numerically for b < bw- Of 
course, the system undergoes a pinning/dcpinning tran- 
sition when crossing the a c — boundary line, but this 
transition is driven by the chemical potential difference 
rather than by the substrate potential and thus it is unre- 
lated to wetting, where phase coexistence of the "liquid" 
and "gas" phases is required (i.e. a = a c ). A very rich 
phenomenology associated with these transitions, have 
been found, however. For b = —4 we find that the non- 
wet phase becomes unstable at a* rj 1.3 and that the 
wet phase becomes unstable at a c —0.064 (as before). 
Consequently, in the range a c < a < a* both phases co- 
exist. This means that if the interface is initially close to 
the wall (n > 1) it remains pinned, while if it is initially 
far from the wall (n ^ 1) it detaches and moves away 
with a constant velocity. Therefore, the system under- 
goes a first-order transition as a function of a. In order 
to establish the phase boundaries we have used the fol- 
lowing criteria. The stability of the pinned phase may 
be characterized by the time r taken by the interface 
to depin in the limit L — > oo. t can be defined as the 
time taken by the last site of the interface to detach, 
/i(x) > or n(x) < 1 Vx. Similarly, we may define 
t as the time characterizing the asymptotic exponential 
decay of (n(t)), where the angular brackets denote av- 
erages of independent runs ( typically 10 5 — 10 6 in our 
simulations). We have verified that both definitions yield 
analogous results. As shown in Fig. for a > a* , r 
saturates with increasing system size and thus the inter- 
face detaches in a finite time. Within the coexistence 
region we have found two different regimes: close to the 
stability threshold of the pinned phase there is a narrow 
stripe 1.22 Js a 1.3 where the detaching time grows 
approximately as a power-law. For a c ^ a ^ 1.22 t 
grows exponentially with L. In both cases r diverges 
as L — > oo, implying that the pinned phase is stable in 
the thermodynamic limit. Due to the very large charac- 
teristic times, we cannot discard the possibility that the 
power-laws are also (asymptotically) exponentials. The 
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study of the asymptotic behavior of the detaching times 
requires longer simulations, beyond our current computer 
capabilities. Finally, the non-monotonic behavior of the 
characteristic times, as well as the step in the curve for 
a = 1.28, may be accounted for by the presence of two dif- 
ferent competing mechanisms as described in jlflj : once a 
site is detached it pulls out its neighbors which, in turn, 
pull out their neighbors in a cascade effect until the whole 
interface is depinned in a time which grows linearly with 
L. This is more likely for small systems, but the proba- 
bility that a site gets detached increases with the system 
size. 

Another way to characterize the power-law regime is 
to analyze the single-site stationary probability density 
function (ss-pdf), defined as the average of n(t) over 
pinned states rather than over all runs. Figure |3| shows 
the unnormalized ss-pdf for different values of a. In the 
exponential-regime (a < 1.22) the histogram exhibits a 
maximum at a pinned state with (n) > 0. In the power- 
law-regime, however, the histogram develops a secondary 
maximum near re = 0, indicating that a fraction of the in- 
terface depins. As a increases, the secondary maximum, 
at zero n, increases while the maximum, at finite re, de- 
creases. At the stability edge (a* « 1.3) the histogram 
changes abruptly into a delta function at n = and the 
pinned phase becomes unstable. 

The differences between the exponential and power-law 
regimes are also observed in a space-time snapshot of a 
numerical solution of ®). In Fig. 0]we plot the stationary 
field re, for a = 1.28, exhibiting patterns characteristic of 
STI [2j| . The main feature of these patterns is the ap- 
pearance of depinned patches (absent for values of a in 
the exponential regime) with a wide range of sizes and 
life-times within the pinned phase. This regime, over- 
looked in previous studies of noncquilibrium depinning 
transitions 0, 0] , seems to correspond to the power- 
law regime described earlier. It is therefore restricted 
to a narrow range between the exponential and the de- 
pinned regimes. This finding is at odds with results of 
previous work claiming that STI is generic in the coexis- 
tence region |22j. 

A typical profile in terms of h is shown in Fig. |SJ The 
depinned intcrfacial regions form triangles with constant 
average slope s. These triangular droplets are similar to 
those described in the discrete model of [3- By taking 
averages of (JSJ), the typical slope, s, of the triangular 
facets is determined through 

\X R \s 2 =a + l 7 (14) 

where Xr is the renormalized non-linear coefficient of the 
KPZ equation. In order to verify equation (|14fl we have 
fixed v = p = D = 1, A = —1 and a = 2.16. Averag- 
ing over 250.000 different triangles, we find an average 
slope s = 1.781, while the value of Xr calculated from 
the tilt-dependent velocity of the depinned interface Q 
yields Xr = —0.9934 from which s — 1.784, in excellent 
agreement with the previously measured value. 

We also studied the size distribution of triangles within 



the power-law regime. Our results correspond to v = p = 
D = 1, A = -1, L = 500, b = -4, and the following val- 
ues of a: 2.12, 2.14, 2.15, and 2.16, and are summarized in 
Fig. a k 2.10 is the boundary between the power-law 
and exponential regimes and the pinned phase is unsta- 
ble for a* w 2.18. The maximum size of the depinned 
regions increases as this instability is approached. Our 
data suggests an exponential dependence on the size of 
the triangular base. This indicates that there is a max- 
imum size for the depinned regions and thus rigorous 
scale invariance (typical of growth driven by a coarsen- 
ing mechanism) of the STI region is ruled out. More ex- 
plicitly, the distribution of triangle sizes, I, is described 
very well by the function exp[3.44(a — 2.176)Z], implying 
that the exponential slopes in Fig. [5] are proportional to 
a — a* . Clearly, triangles with a base less than ^2,3 
cannot be visualized due to the discretization of equa- 
tion ©. A simple extrapolation indicates that the tri- 
angles become imperceptibly small for a rs 2.05, in good 
agreement with the value obtained for the boundary be- 
tween the power-law and exponential regimes. Therefore, 
we cannot rule out the possibility that the triangles are 
ubiquitous throughout the coexistence region (although 
not always visible in a discrete numerical simulation) in 
which case the power-law dettaching times should turn 
into exponentials for large enough times and system sizes. 
In this case the force exerted by the non-linear KPZ term 
on the triangular facets against the direction of growth, 
at late times, guarantees the stability of the pinned phase 
fl8| throughout the finite coexistence region. 

Finally, we study the phase-coexistence regime in a 
two-dimensional system to check whether the triangular 
patterns survive in higher dimensionalities. In particu- 
lar, we consider a system size 100 x 100 and take a = 2, 
b = —4, well within the coexistence region. We find 
structures as those shown in Fig. [7| the triangles becom- 
ing pyramids. Note that the edges of the pyramid bases 
are parallel to the axes of the discretization-lattice. This 
suggests that the pyramids are lattice artefacts and that 
a continuum system may exhibit conical structures. 



V. CONCLUSIONS 

We have investigated a continuum model for nonequi- 
librium wetting transitions. The model consists of a KPZ 
equation in the presence of a short-ranged substrate po- 
tential, and is the most natural non-equilibrium exten- 
sion of the interface displacement models used in equi- 
librium wetting. It can be mapped into a multiplicative 
noise problem, enabling simple theoretical calculations at 
the mean-field level. Numerical simulations reproduce a 
phase diagram analogous to that obtained within mean- 
field, including first as well as second-order phase tran- 
sitions. In particular, we have found complete wetting 
and critical wetting transitions, as well as a finite area in 
the temperature-chemical potential phase diagram where 
pinned and depinned phases coexist. This finite coexis- 
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tence region allows us to define critical wetting along a 
range of paths that are, however, characterized by the 
same critical exponents. Within this area we identified 
two regimes. In the first, the lifetime of the pinned phase 
grows exponentially with increasing system size and its 
ss-pdf is bell-shaped. The second one exhibits STI, life- 
times consistent with a power-law, and a double-peaked 
ss-pdf. The main feature of the latter regime is the pres- 
ence of triangular structures that have been characterized 
by their slopes and size distributions. 

An interesting open problem is that of the equilibrium 
limit of non-equilibrium wetting. The Cole-Hopf trans- 
form precludes the limit A = to be studied using this 
method. Moreover, we have noted how the behavior of 
the EW equation in the presence of a wetting potential 
differs from the weak-noise regime of the MN equation. 
This leaves the crossover to equilibrium wetting an open 
challenge. In addition, the effects of long-ranged poten- 
tials on the phenomenology described here remain to be 



investigated. 

Finally it would be extremely interesting to develop 
experiments in order to explore the rich, non-equilibrium 
phenomenology described in the previous sections; liquid- 
crystals are in our opinion good candidates for this. It 
is our hope that this work will stimulate experimental 
studies in this direction. 
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FIG. 1: Mean-field phase diagram and typical solutions of 
the equation F((n}) — (n) (temperature in units such that 
k B = 1. 
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FIG. 2: Top: linear correlation coefficient for Inn as function 
of In 1 6 — bw\ for different values of bw- The maximum gives 
the best estimate of bw = —0.70 ±0.01. Bottom: log-log plot 
of Inn as function of In |6 — bw\ in the vicinity of the critical 
point. The line is a least-squares fit; from its slope the critical 
wetting exponent is found to be 1.20 ± 0.02 
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FIG. 3: (A) Characteristic depinning times and (B) ss-pdf for 
various representative values of a. 




FIG. 4: Configuration in the n-representation for a = 1.18 
and b = —4. Depinned regions (n < 1) are colored in dark- 
grey and pinned ones (n > 1) in light-gray. 1000 time slices 
are depicted at intervals of 50 time uints. The system size is 
L = 500. 
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FIG. 5: Instantaneous configuration of the interface for time 
slice 400 (marked with a line in Fig. . Parameters as in Fig. 




FIG. 6: Distribution of triangles as a function of the size of 
the triangular base, for a = 2.12, 2.14, 2.15 and 2.16. 



FIG. 7: Snapshot of an interface configuration for a 100 x 100 
system (not all the substrate is shown) and parameters a = 2 
and b — —4. 



